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ABSTRACT 

Ultra compact dwarf galaxies (UCDs) have similar properties as massive globular 
clusters or the nuclei of nucleated galaxies. Recent observations suggesting a high 
dark matter content and a steep spatial distribution within groups and clusters provide 
new clues as to their origins. We perform high-resolution iV-body / smoothed particle 
hydrodynamics simulations designed to elucidate two possible formation mechanisms 
for these systems: the merging of globular clusters in the centre of a dark matter 
halo, or the massively stripped remnant of a nucleated galaxy. Both models produce 
density profiles as well as the half light radii that can fit the observational constraints. 
However, we show that the first scenario results to UCDs that are underluminous 
and contain no dark matter. This is because the sinking process ejects most of the 
dark matter particles from the halo centre. Stripped nuclei give a more promising 
explanation, especially if the nuclei form via the sinking of gas, tunneled down inner 
galactic bars, since this process enhances the central dark matter content. Even when 
the entire disk is tidally stripped away, the nucleus stays intact and can remain dark 
matter dominated even after severe stripping. Total galaxy disruption beyond the 
nuclei only occurs on certain orbits and depends on the amount of dissipation during 
nuclei formation. By comparing the total disruption of CDM subhaloes in a cluster 
potential we demonstrate that this model also leads to the observed spatial distribution 
of UCDs which can be tested in more detail with larger data sets. 
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1 INTRODUCTION 

A new population of subluminous and extremely compact 
objects have been recently discovered in cluster and group 
environments (Drinkwater et al. 2000; Phillips et al. 2001) 
These ultra compact dwarf galaxies (hereafter UCDs) are 
dynamically distinct systems having intrinsic sizes ~ 100 pc 
and absolute magnitudes in the _B-band in the range from 
— 13 to —11 placing them in the lower range of dwarf galaxy 
luminosities (Mateo 1998). A number of different scenarios 
have been proposed for origin of UCDs (e.g. Evstigneeva 
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et al. 2007 and references within): (a) They are simply very 
big and luminous globular clusters, (b) they are nucleated 
dwarf galaxies, (c) they are the resulting objects of the co- 
alescence of several globular clusters, (d) they are the rem- 
nants of stripped disk galaxies. In this work we compare the 
last two scenarios. 

Oh and Lin (2000), Oh, Lin and Richer (2000) as well as 
Fellhauer and Kroupa (2002) investigated the third forma- 
tion scenario. The latter authors performed A'^-body simula- 
tions of the merging of several globular clusters and argued 
that the resulting object resembled a UCD. Their simula- 
tions initiated with a supercluster, an accumulation of up 
to fifty globular clusters, orbiting inside a host galaxy on 
fairly elliptical orbits and with apo-galactic distances around 
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20kpc. The globular clusters merged within the superclus- 
ter giving rise to an object with values for surface brightness 
and absolute bolometric luminosity comparable to UCDs. 

The last formation scenario has been discussed in Oh, 
Lin & Aarseth (1995) as well as Bekki et al. (2001, 2003). 
The latter authors performed numerical simulations of the 
dynamical evolution of nucleated dwarf galaxies orbiting in- 
side NGC 1399 and the Fornax cluster. Adopting a plausible 
scaling relation for dwarf galaxies, Bekki et al. found that 
the outer stellar envelope of a nucleated dwarf was totally re- 
moved by tidally stripping over the course of several passages 
from the central region of their host. The nucleus is so dense 
that it will always survive the tidal field of a group or clus- 
ter potential. By construction in the initial conditions, the 
size and luminosity of the remnant were similar to those ob- 
served for UCDs and the host galaxy was a Plummer model 
halo which has a constant density core and therefore easy 
to tidally disrupt leaving behind the central nucleus which 
would be associated with a UCD galaxy. 

The main goal of the present study is to shed light into 
the formation of UCDs investigating the last two formation 
scenarios in more detail and with more realistic initial con- 
ditions. First, we adopt the Fellhauer and Kroupa (2002) 
model and combine it with the cold dark matter model 
(CDM) paradigm. We assume that globular clusters form 
and subsequently orbit around dark matter haloes having 
masses comparable to that of the Fornax dwarf spheroidal 
and containing no other baryonic matter. Due to dynam- 
ical friction, these globular cluster would spiral in towards 
the centre of the halo. We perform coUisionless A''-body sim- 
ulations of this process and show that the object resulting 
from the coalescence of the globular clusters at the halo cen- 
tre resembles a UCD. As we demonstrate below, the above 
model suffers from two major drawbacks which rule out its 
applicability for UCD formation. 

The first issue is related to the observed high M/L ratio 
(between 6 and 9 in solar units) recently reported for UCDs 
(Ha§egan et al. 2005). One has to note here that the M/L 
ratio of UCDs is still debated (Micheal Drinkwater, private 
communication, compare Evstigneeva et al. 2007 who quote 
a M/L for UCDs 3 and 5 in solar units). Such high values can 
probably not be achieved by the above mechanism. This is 
because sinking globular clusters will expel most of the dark 
matter particles from the halo centre (El-Zant, Shlosman 
& Hoffman 2001; Merritt et al. 2004; Goerdt et al. 2006). 
If there is no dark matter in UCDs, Fellhauer and Kroupa 
(2006) describe a possible way to enhance the mass-to-light 
ratios of UCDs through tidal interactions. The second prob- 
lem with this scenario is related to the total luminosity of 
an UCD. At least today, dark matter halos with a sufficient 
number of globular clusters to produce such bright objects 
are very rare (Sharina, Puzia & Makarov 1996; Boker et al. 
2002). 

Second, we examine the scenario of Bekki et al. (2001, 
2003) which is based on the hypothesis that UCDs are rem- 
nants of stripped nucleated galaxies. This mechanism has 
also been discussed by Kazantzidis et al. (2003). We test 
this model using SPH simulations of a low mass galaxy 
which forms a nucleus via gas inflow to the inner ~ 100 
parsecs. Once this galaxy is placed on a critical disrupting 
orbit within a cluster potential, the surviving nucleus is in 
excellent agreement with the latest observational constraints 



for UCDs, including their dark matter content and spatial 
distribution. 

The paper is organised as follows. In section 2 we de- 
scribe the globular cluster numerical simulations and com- 
pare the properties of the resulting object with those of 
UCDs. Section 3 contains results from the simulations of 
tidal stripping of disk galaxies inside a host cluster potential. 
Finally, in Section 4 we summarise our main conclusions. 



2 THE GLOBULAR CLUSTER SIMULATIONS 

The globular cluster simulations were performed with Pkd- 
GRAv2 a multi stepping, parallel A'"-body tree code (Stadel 
2001). We create an NFW (Navarro, Frenk & White 1996) 
halo, employing the technique developed by Kazantzidis, 
Magorrian & Moore (2004), which has the following den- 
sity profile: 



p(r) 



Po 



r/rs [1 + (r/rs)] 



(1) 



In our case po ~ 242 Mo/pc"' and rs = 1.5 kpc. The halo has 
a virial mass of 1.5 x 10^ Mq . The concentration parameter is 
20 which is the typical value for halos of this initial mass. To 
increase mass resolution in the region of interest, we divide 
the halo into three shells (Zemp et al. 2007) each of which 
contains 10"" particles. The innermost shell has 100 pc radius. 
The second shell is between 100 and 500 pc while the third 
shell contains the rest of the halo. The softening lengths 
for these shells are 1, 10 and 100 pc respectively. This shell 
model allows us to resolve the detailed kinematics within 
the central few parsecs whilst retaining the global structure 
of the halo out to its virial radius of 29.39 kpc. 

We use ten globular clusters consisting of 10^ particles 
and being represented by the King model (King 1966; Michie 
1963; Michie & Bodenheimer 1963) 



(2) 



p(*) = pi exp — erf 



/ 4* r 2* 
P^\ — 2 1 + TT 



Each globular cluster is constructed with a Wo = 'I'(0)/(7^ 
parameter of 6, a total mass of 4.2 x 10' Mq, a central ve- 
locity dispersion of llkm/s, and an absolute magnitudes 
of -8.5, assuming a mass to light ratio of 2 (Boker et al. 
2004; Walcher et al. 2005) . We use 1 pc for its gravitational 
softening length. The ten globular clusters are initially dis- 
tributed within the halo between 20 pc and 1000 pc. They 
are spatially distributed according to p{r) oc r~*, which is 
in agreement with observations of globular cluster popula- 
tions in dwarf galaxies (Sharina, Puzia & Makarov 1996). 
We randomly place the globular clusters on circular orbits 
(We also tried different orbital and spatial distributions in- 
cluding r~^ and r" with similar results as in the standard 
r case.). 

The dwarf halo which contained the globular clusters 
was placed on a static NFW potential which corresponds 
to a cluster halo with virial mass of Mvir = 10^'* M© and 
concentration of c = 6. The halo was placed on an initial 
distance of 100 kpc from the centre of the potential and was 
given a velocity of 650 km/s, respectively. These values result 
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Figure 1. The diflferent density profiles of tlic UCD and the 
corresponding darlc matter distributions, after the ten globular 
clusters merged, but before and after it has been put into an 
external potential. On top of the density profile of the final stellar 
distribution we show the best fit plummer model. 



to an eccentric orbit with an apocentre of 500 kpc. For a 
second simulation we put fifty of these globular cluster in 
the aforementioned halo and let them spiral into the centre. 
This second simulation identical to the first one, just with 
fifty instead of ten globular clusters. 

All ten globular clusters in our first simulation merge 
to a single object within 0.4 Gyr. The associated timescale 
of this evolution can be calculated using the Chandrasekhar 
dynamical friction formula (Chandrasekhar 1943). 



dr 
df 



47rlnA(r)p(r)G^MGcr 
Vc{r)d [rvc{r)] /dr 



erf 



\/2a(r; 



2v^{r) 



27r cr(r) 



exp 



-vl{r) 



2cr2(r) 



(3) 



where vdr) is the circular speed at radius r, In A(r) is the 
Coulomb logarithm, which we assume to be 4.0, Mac is the 
mass of one globular cluster. p{r) is the density of the dark 
matter halo at radius r according to equation (1) and (j{r) is 
the one-dimensional velocity dispersion of the halo. Solving 
this equation numerically gives thic ~ 0.4 Gyr, which is very 
similar to the value we directly get from the simulation. 

The central compact object which results from the co- 
alescence of the globular clusters has a three-dimensional 
density profile which is shown in Fig. 1. The profile is very 
similar to ones derived from observations (De Propris et al. 
2005) as well as previous simulations (cf. Fig. 5 in Fellhauer 
& Kroupa (2002)). The density profile of the baryonic mat- 
ter only (i.e. just the material, which has been in the original 
globular clusters) can be well-fitted with a Plummer model 
(Plummer 1911): 

5/2 



p(r) = 



3M 






(4) 



In our case M = 1.45 x lO'' Mq and b = 7.26 pc. The 
resulting object has a half mass radius of 25 pc and a central 



Figure 2. Projected surface density profile of the UCD, after the 
ten globular clusters merged, but before it has been put into an 
external potential. The profile of the luminous matter together 
with the best fit plummer model is plotted. 



velocity dispersion of 20 km/s. The projected surface density 
profile can be seen in Fig. 2. The fit is good in the inner parts 
and reasonable in the outer parts. The equation for the fit 
has been found by Evans & An (2005): 

A 



S(r) 



7r(l -FrVc2)' 



(5) 



Our fit gives A = 3.3 x lO^Mo/pc^ and c = 8.7pc. It 
looks fairly similar to observations and previous simulations 
[cf. Fig. 5 in Ha§egan et al. (2005) or Fig. 6 in Fellhauer 
& Kroupa (2002)]. If we assume a mass to hght ratio of 
three, we expect an absolute visual magnitude of -11 for the 
final object. For comparison, a typical UCDs found in For- 
nax has a central velocity dispersion of 22 km/s, an absolute 
magnitude of -12 and a half mass radius of 25 pc (Ha§egan et 
al. 2005; Drinkwater et al. 2003). In order to form this sys- 
tem by merging globular clusters one would need to consider 
many more globulars to reach the observed luminosity. This 
turns out to be a crucial point of this formation scenario 
because the presence of more than ten globulars in such a 
halo is extremely rare (Sharina, Puzia & Makarov 1996). 

In our second simulation, the fifty globular clusters 
merged completely within a few million years. The result- 
ing object had a central velocity dispersion of 37.3 km/s and 
an absolute magnitude of -12.75, assuming a mass to light 
ratio of three. Its position is marked in Fig. 3 and roughly 
matches a UCD. 

Interestingly, all dark matter particles are expelled from 
the core of the final object and the inner initial cusp has 
been turned into a nearly constant density core. Indeed, in 
our ten globular cluster simulation not even a single dark 
matter particle is found within the inner 5pc. This effect 
has been discussed in more detail in Read et al. (2006) 
and constitutes the second reason for excluding the merg- 
ing of globular clusters as a possible formation mechanism 
for UCDs. In other words, the observed M/L ratio cannot 
possibly be reached without dark matter. Only if the initial 
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Figure 3. Central velocity dispersion versus absolute visual mag- 
nitude for various stellar systems. Observations are plotted to- 
gether with the results from our simulations. Data taken from 
Ha§egan et al. (2005); Evstigneeva et al. (2007) and references 
therein. 



globular clusters contained cuspy dark matter distributions 
would such a scenario stand a chance of working, however 
current observations do not support this idea in the galactic 
globular clusters. 

Keeping in mind the observational uncertainty of the 
M/L ratios of UCDs, the mechanism proposed by Fellhauer 
& Kroupa (2006) to enhance the M/L ratio of UCDs as 
well as resolution limits and other shortcomings of properly 
determining the M/L ratio of our simulations in the next 
section, the reader should note that the major issue and 
stronger argument against this formation scenario is the to- 
tal central luminosity. 



3 ARE UCDS STRIPPED DISK GALAXIES? 

The second scenario we consider for the formation of UCDs 
is that they are the remnants of stripped disk galaxies. For 
the numerical experiments performed in this section, we 
build galaxy models using the following procedure. First, 
we set up a spherical equilibrium NFW halo with struc- 
tural parameters consistent with predictions of the standard 
ACDM model (Kazantzidis, Magorrian & Moore 2004) . We 
include a gaseous component of mass equal to a fraction /b 
of the total halo mass. The gas has originally the same ra- 
dial distribution and a temperature profile such that it is 
in hydrostatic equilibrium for an adiabatic equation of state 
(EOS). 

The gas component has a specific angular momentum 
distribution and spin parameter consistent with values found 
for dark matter haloes within cosmological N-body simu- 
lations (Bullock et al. 2001; Maccio et al. 2007). We con- 
structed a dark plus gaseous halo model with parameters 
that are expected to produce disks similar to the Local 
Group galaxy M33 which has a small dense stellar nucleus. 
The model parameters were; Mvir = 5 x 10^^ Mq, r^ir = 167 



kpc, Wvir = 115 km/s, c — 6.2, baryonic fraction /b = 6% 
and spin parameter A = 0.1. 

The value of the concentration may seem a little bit low, 
when comparing it to the results of Maccio et al. (2007) but 
it is the best fit for the M33 galaxy (Corbelli 2003) which is 
the possible prototype progenitor for UCDs, we want to test 
here. The value of the spin parameter is what one gets from 
cosmological dark matter simulations. This net spin results 
from adding many large, almost randomly oriented angular 
momenta of individual dark matter particles. Uniformly ro- 
tating gas halos are normally not observed in cosmological 
simulations, gas is accreted pretty much like the dark mat- 
ter, leading to random, large bulk motions and turbulence 
(Wise & Abel 2007). So one has to note here that this setup 
is designed to produce cored disks for the subsequent exper- 
iments and not necessarily an attempt to follow disk and 
core formation realistically. 

The hot gaseous halo is resolved with 2 x 10® particles 
of equal mass ~ 2 x 10* M0 . We sampled the dark matter 
halo with 2.2 x 10® particles having variable masses with 
the resolution increasing towards the centre of the system 
(Zemp et al. 2007). With a single-particle model one would 
need about ten million particles to reach a comparable res- 
olution in the central region. This allows us to simulate the 
central dark matter cusp with softening of 100 pc (this is 
the same for all dark matter and gas particles) and particle 
mass of ~ 4.4 x lO** Mq. The detailed description of the 
initial conditions and results of the evolution of the disk are 
presented in Kaufmann et al. (2006) and Kaufmann et al. 
(2007). 

We use the parallel Tree-j-SPH code Gasoline (Wads- 
ley, Stadel & Quinn 2004) , which is an extension of the pure 
N-Body gravity code Pkdgrav developed by Stadel (2001). 
It uses a 32 particles smoothing kernel and an artificial vis- 
cosity which is the shear reduced version (Balsara 1995) of 
the standard Monaghan (1992) artificial viscosity. Gasoline 
uses a spline kernel with compact support for the softening of 
the gravitational and SPH quantities. The energy equation is 
solved using the asymmetric formulation, which is shown to 
yield very similar results compared to the entropy conserving 
formulation but conserves energy better (Wadsley, Stadel & 
Quinn 2004). The code includes radiative cooling for a pri- 
mordial mixture of helium and (atomic) hydrogen. Because 
of the lack of molecular cooling and metals, the efficiency 
of our cooling functions drops rapidly below 10 000 K, but 
we adopt a temperature fioor of Tj = 15 000 K, to crudely 
mimic the effect of heating sources such as supernovae ex- 
plosions and radiation (e.g. ultraviolet) backgrounds. These 
simulations were performed without star formation. We do 
not expect them to suffer from the SPH problems pointed 
out by Agertz et al. (2007) because for the disk forming 
part of the simulations we are adopting a large temperature 
fioor and do not aim to simulate the actual multiphase ISM, 
where the problems would occur. The gas cools slowly to 
form a disk and stars are formed at a rather low thresh- 
old, therefore there is problem with our simulations (Oscar 
Agertz, private communication). For a detailed discussion on 
the performance of SPH for the stripping part of the simu- 
lations we defer the interested reader to McCarthy (2007). 

Due to cooling the hot gas halo loses its hydrostatic 
equilibrium quite quickly. An inner gaseous disk rapidly 
forms out of cooling gas coming from a nearly spherical re- 
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Figure 4. The logarithmic surface density of the M33 gas dislc is 
plotted after 2.5 Gyr. The nucleus in the centre is clearly visible. 
The dark matter distribution is plotted on top for comparison. 



Figure 5. The dark matter distributions of the disk model in 
various stages of its evolution, together with the final baryonic 
(stellar) mass distribution. 



gion close to the halo centre (within ~ 10 kpc). After 2.5 
Gyr of evolution the disk attains a near exponential sur- 
face density profile over a large fraction of its extent, except 
within a few hundred parsecs from the centre, where gas in- 
flows along an inner bar produce a dense nucleus which can 
be seen in Fig. 4. 

The nucleus has a mass of --^ 3 x 10^ Mq, which is 
~ 0.6% of the total disk. Central dense stellar nuclei have 
been seen in several late-type spiral galaxies, for example, in 
M33 (Regan & Vogel 1994), but also in other late-type spi- 
rals (CaroUo 1999; CaroUo, Stiavelli & Mack 1998). Recent 
observations by Ferrarese et al. (2006) show that they are 
observed in 50 - 80% of low- and intermediate-luminosity 
galaxies. 

We put the galaxy into the same artificial NFW po- 
tential resembling the Virgo cluster that we used before. We 
start with the disk at pericentre, 100 kpc far from the centre 
of the static potential, and give it a perpendicular velocity 
of 260 km/s, such that it should reach an apocentre at 30 
kpc. The plane of the disk is tilted by 45° with respect to 
the plane of the orbit. On this orbit the outer galactic halo 
and gas are quickly stripped from the nucleus. 

The resolution of our hydrodynamic simulation is not 
sufficient to allow a meaningful comparison with the obser- 
vations. However the half light radius, luminosity and dark 
matter content all agree well with observations up to the 
given resolution limits (see Fig. 5), the mass to light ratio 
within the optical radius (~ 50pc) is about 4. and the UCD 
is dark matter dominated up to the resolution limit. 



3.1 Defining the critical disruption zone and tiie 
affects of baryons on survival versus 
disruption 

If the orbit of the galaxy does not penetrate the centre of the 
cluster then the disk would not be stripped and the object 
would not be classified as a UCD. If the progenitors are 



nucleated disks which are tidally transformed into nucleated 
spheroidals dE(N) via galaxy harassment (Moore et al. 1996) 
then we might expect a transition region from central UCDs 
to outer dE(N). Since UCDs are all located close to the 
centres of clusters and groups this allows us to explore the 
conditions under which complete disruption occurs and how 
that is affected by the amount of dissipation (steepening of 
the central potential). 

In what follows, we performed an extensive set of sim- 
ulations to identify the orbits on which the disk galaxy is 
completely stripped and on which it remains intact. The 
first case would correspond to a UCD and the second to a 
dE(N). The central potential of the SPH galaxy is deepened 
significantly due to dissipation. The object will be harder to 
completely disrupt as compared to the initial dark matter 
halo of the model for example. We therefore make a sys- 
tematic study of the orbits that lead to disruption versus 
survival for the nucleated M33 galaxy model and for the ini- 
tial uncompressed dark matter halo - the latter simulations 
define the optimum scenario for complete disruption. 

These simulations all start at their respective apocen- 
tres and run for the same physical time (5 Gyr) in the cluster 
potential, and explore different orbital apocentres and peri- 
ccntres. If stripping of disk galaxies is the correct mechanism 
for forming UCDs and dE(N), there must be a very sudden 
transition between these two cases. This is because we do 
not observe an intermediate object with features lying in 
between these very distinct two types of galaxies. 

Fig. 6 presents surface brightness profiles of the disk 
galaxy after orbiting the static cluster potential for 5 Gyr 
on various orbits. A sudden transition is clearly seen in this 
plot: almost two orders of magnitude in surface brightness 
lie between the (150 - 30) and the (125 - 25) orbit at the 
radii of interest, so the disk either survives or is completely 
disrupted. 

Fig. 7 presents a scatter plot of apocentres versus peri- 
centres for all cases we simulated along with indications 
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Figure 6. Surface brightness profiles of the stripped dislc after or- 
biting an artificial Virgo cluster potential for 5 Gyr. The numbers 
next to the hyphen in the legend give the apo- and pericentres 
of the various orbits in kpc and the numbers within the brack- 
ets give the number of orbits which were completed during the 5 
Gyrs. The solid line indicates the dividing line between survival 
and disruption. 



about the fate of the orbiting system. The dividing Une be- 
tween disrupted and surviving disks is given by 



^rapo + B 



(6) 



with A = -0.32 and B = 68. 

We also repeated the above stripping experiments with 
the dark matter halo only without the gas disk. Results of 
these experiments are included in Fig. 7. In this case the 
transition is not as sudden as in the disk case therefore we 
define a halo to be completely disrupted, if less than 0.1 % of 
the virial mass is still bound. Defining complete disruption 
by the bound mass getting below a certain threshold is the 
most objective way to do so. The value of 0.1%, may seem 
arbitrary, but it is actually the value, which gave the most 
convincing values for the given resolution and algorithm. To 
estimate the mass which is still bound we used the group 
finding algorithm Skid (Stadel 2001). There is a clear sep- 
aration between survival and disruption that depends only 
on the apocentre and pericenter. The theoretical explana- 
tion for the orbits that lead to complete disruption is fairly 
intricate and will be discussed in a later paper. 

The dividing line between disrupted and surviving 
haloes is given roughly by two straight lines: 



Ar^po + B 
Crapo + D 



(7) 



with A = -0.11, B = 84.0, C = -0.8, D = 335.0, with the 
transition between them at r = (_D — _B)/(A— C) ^ 370 kpc. 
Given these empirically determined survival/disruption 
curves, we can sample orbits from a cosmological cluster 
mass CDM halo and determine the expected radial profile 
of completely disrupted haloes (UCDs) in each of these two 
cases. From the cosmological simulation we chose the one 
cluster which was not perturbed by neighbouring clusters. 



Figure 7. Disrupted and non disrupted gas disks and dark matter 
halo nuclei as a function of apo- and pericentre for all the disk 
and the dark matter only simulations we have done. The solid 
line is the bisector and the two dotted ones are the dividing lines 
between disruption and survival. 



had only a single nucleus and came closest in Mvir and rvir 
to the Virgo cluster. 

We evolved the halo for 5 Gyr and binned all parti- 
cles according to their projected distance from the centre at 
the last simulation output. For each bin we determined the 
fraction of particles whose orbits during the 5 Gyr of evo- 
lution lies below the respective dividing line in Fig. 7. The 
orbits have been determined in a way that the pericentre 
is the closest point the particle gets to on all snapshot and 
the apocentre is the furthest point. In Fig. 8 this fraction is 
then compared to its observed counterpart in nature. This 
is the fraction of the number of UCDs to the number of pro- 
genitors, which is the sum of the number of UCDs and the 
number of dE(N)s, so: 



/uCD 



nvcn 



nucD 



T-progcnitor l-UCD + fldE(N)s 



(8) 



Results from the pure dark matter stripping experi- 
ments correspond to the upper cosmological prediction (The 
final radial distribution of UCDs is more extended since the 
dark matter halo is not compressed due to dissipation and 
disrupts at larger clustercentric radii) . The SPH disk galaxy 
stripping simulations correspond to the lower cosmological 
prediction. We have to remind the reader again that our pro- 
totype progenitor was M33 that has a concentration, which 
might look a bit low according to Maccio et al. (2007). 
Therefore the subhalo disruption and the UCD fractions 
might be slightly overestimated. However, one could image 
the progenitor being heavier, therefore less concentrated and 
thus showing the same disruption behaviour as the haloes we 
present. Unfortunately we did not have the computational 
resources to explore parameter space any further. The obser- 
vational data have been obtained in the following way. The 
projected distances of UCDs and dE(N)s in arcsec (Jones 
et al. 2006) are converted into kpc assuming that the Virgo 
and Fornax cluster lie at distances of 16 and 20 Mpc, re- 



© 2007 RAS, MNRAS 000, 1-7 



Formation of UCDs and dE(N)s 7 



Q 
O 

3 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




upper cosmological prediction 

lower cosmological prediction 

observational data 




50 100 150 200 250 

r/kpc 



Figure 8. The ratio of the number of UCDs to the number of 
progenitors, which is sum of the number of UCDs and the number 
of dE(N)s, as a function of radius. The upper cosmological predic- 
tion comes from the disruption simulations using the pure dark 
matter subhalo whilst the lower cosmological prediction comes 
from the SPH simulations with the dissipated gas disk. Observa- 
tional data is taken from Jones et al. (2006). 



spectively. In a ACDM universe, haloes are approximately 
self similar so we can compare results by scaling everything 
with respect to the virial radius. After rescaling, we binned 
the occurrence of the UCDs and dE(N)s according to their 
rescaled projected distance to their mutual galaxy cluster 
centre into common bins. The content of each bin is associ- 
ated with a Poissonian error. Finally, we divide the number 
of UCDs in each bin by the sum of the numbers of UCDs and 
dE(N)s and propagate the errors accordingly. From compar- 
ing our two cosmological prediction with the data, we con- 
clude that (a) both UCDs and dENs originate from the same 
progenitors and (b) this progenitor must have properties in 
between the disk galaxy and the pure dark matter halo we 
adopted in our simulations. One must bear in mind though 
that we only have position data for 15 UCDs from the lit- 
erature. Nonetheless, and limitations notwithstanding, our 
cosmological predictions agree well with the data but further 
observations will allow a better test of this model. 



4 CONCLUSIONS 

UCDs have similar properties as massive globular clusters 
or the nuclei of nucleated galaxies. Recent observations of 
a high dark matter content and their steep spatial distri- 
bution within groups and clusters give us new clues as to 
their origins. We perform A^-body simulations and compare 
two possible mechanisms for their formation: the merging 
of globular clusters in the centre of a dark matter halo, or 
the massively stripped remnant of a nucleated galaxy. Our 
simulations reveal, that a swarm of ten as well as one of 
fifty globular clusters born in a Fornax dwarf spheroidal size 
dark matter halo will normally evolve to an object, which 
has the same density profile as an UCD. We performed a 



second series of simulations in which a disk galaxy experi- 
ences mass loss processes inside a cluster environment. The 
disk is entirely stripped and the remaining nucleus exhibits 
all of the observed UCD properties. Both models produce 
density profiles as well as the half hght radii that can fit the 
observational constraints very well. 

However, we show that the first scenario produces UCDs 
that are underluminous and contain no dark matter - the 
sinking process ejects most of the dark matter particles from 
the halo centre. The stars, which consist the globulars, would 
replace the dark matter particles and expel them from the 
centre of the halo (El-Zant, Shlosman & Hoffman 2001; Mer- 
ritt et al. 2004) lowering the M/L ratio beyond observational 
constraints (Ha§egan et al. 2005). The other drawbacks of 
this model is that it is very unlikely to have so many glob- 
ular cluster in a halo of that size. It becomes apparent in 
Fig. 3 that it is not sufficient to merge only ten globulars: 
approximately fifty are necessary. However such a halo may 
not exist (Sharina, Puzia & Makarov 1996). 

Stripped nuclei give a more promising explanation, es- 
pecially if the nuclei form via the sinking of gas, funneled 
down inner galactic bars, since this process enhances the 
central dark matter content. Even when the entire disk is 
stripped away, the nucleus remains intact and can be dark 
matter dominated. The total disruption of the galaxy be- 
yond the nuclei only occurs on certain orbits and depends on 
the amount of dissipation during nuclei formation. By com- 
paring the total disruption of CDM subhaloes in a cluster 
potential we show that this model also leads to the observed 
spatial distribution of UCDs. This scenario can be tested in 
more detail with larger data sets. 
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